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ABSTRACT 

The nearby galaxy M33 hosts the most luminous steady X-ray source in the Local Group. 
The high spatial resolution of Chandra allows us to confirm that this ultra-luminous X-ray 
source is within the nucleus and rule out at the 4.6a level a previously proposed possible 
counterpart located 1" away. The X-ray spectrum is well fitted by a disc blackbody with 
fcX- n =1.18±0.02keV and i? in (cosf?)°- 5 =57±3 km, consistent with earlier results from ASCA 
and BeppoSAX. The source flux is steady between 1 and 10 4 s (<3% rms variability, 0.5- 
10 keV) with a 0.5-10 keV luminosity of 1.5 x 10 39 erg s _1 . The X-ray properties and an 
association with a radio source are reminiscent of the galactic microquasar GRS 1915+105. 

Key words: galaxies: individual (M33) — galaxies: nuclei — Local Group — X-rays: galax- 
ies. 



1 INTRODUCTION 



The nearby galaxy M33 (d « 795 kpc; van den Bergh 1991) hosts 



the most luminous steady X-ray source in the Local Group. M33 X- 
8 dominates the X-ray emission from the galaxy with a 1- 10 keV lu- 
minosity of about 1.2-10 39 erg s" 1 , ten times more than that of the 
other X-ray sources in the field. The ultra-luminous X-ray source 
(ULX) has been detecte d at this level since its d iscovery by Ein- 
stein . (^ong et al. 198~j]; |vlarkert & Rallis 1983|; Gottwald et al. 



1 987; [Peres et al. 1989fc [Takano et al. 1994j |Long et al. 1996j; P ar- 
mar et al. 2001). Observations obtained with the ROSAT High Res- 
olution Imager (HRI) showed that X-8 is point-like and coincident 
with the optical n ucleus of M33 at the HRI 5" resolution ( Schulman 
& Bregman 1995). These also reveal ed a 106 day ^20 % modula- 
tion of the X-ray flux from M33 X-8 ( frjubus et al. 1997| ). 

The nature of the X-ray emitting object is rather puzzling. The 
measured velocity dispersion of the nucleus places an upper limit 
of 1500 Mq on the m ass of a putative massive central black hole 
( |Gebhardt et al. 2001 ). A low luminosity active galactic nucleus 
(AGN) is therefore unlikely. Indeed, the UV/optical spectrum of the 
nucleus shows no signs of AGN-type activity and is best modelled 
by a combination of t wo starbursts at 40 Myrs and 1 Gyrs (L ong, 
Charles & Dubus 2002 and references therein). 

The 106 day modulation suggests a single object is respon- 
sible for most of the X-ray emission, perhaps a persistent X-ray 
binary with a 10 Mq black hole radiating at the Eddington luminos- 
ity. Hubble Space Telescope observations place a limit of 18% on 
the contr ibution of a point so urce to the total far-UV flux from the 
nucleus (Dubus et al. 1999b). The implied ratio L opt /Lx ~ 0.05 



would be in line with expectations from a low-mass X-ray binary. 
The counterpart would then be mostly lost in the glare of the nu- 
cleus. However, there remains a distinct possibility that M33 X-8 
is associated^vvitlijiUVJjngJW^s^an^ 

et al. 1999b). The star has a 09III spectrum with no outstanding 



Table 1. 


Chandra Observations of M33 






Det. 


Obsid 


Seq.# 


Date 


Exp. 


offset 


ACIS-I 


1730 


600145 


7/12/2000 


49.4 ks 


0' 


ACIS-S 


786 


600089 


8/30/2000 


46.2 ks 


0' 


ACIS-S 


787 


600090 


1/1 1/2000 


9.2 ks 


7.7' 



Det. is the detector used; exp. is the total integration time; offset is the dis- 
tance of M33 X-8 from the telescope axis. The CCD frame time was 3.2 s 
for the first two observations and 1.0 s for the last, offset observation. 



feature linking it to the powerful X-ray source but cannot be ru led 
out at the ROSAT HRI resolution ( |Long, Charles & Dubus 2002| ). 

Most studies have assumed M33 X-8 is associated with the 
nucleus and used this to register the X-ray images. Chandra offers 
the opportunity to locate M33 X-8 to unprecedented accuracy and 
settle whether this source is within the nucleus or not. Three obser- 
vations of M33 are available in the Chandra archive (Table 1). We 
used the first two to accurately position M33 X-8 with respect to 
HST Planetary Camera (PC) images of the nucleus (§2). The third, 
non piled-up observation is used to study the X-ray spectrum and 
timing behaviour of M33 X-8 (§3 and §4). Our results are discussed 
in the last section. 



2 ASTROMETRY 

The 1-2" precision of the Chandra pipeline aspect solution^] is in- 
sufficient for our purposes. To improve on this, we identified X-ray 
sources within the field-of-view and cross-correlated these to opti- 
cal, IR and radio sources with accurately known positions. 



http://cxc.harvard.edu/cal/ASPECT/celmon/index.html 
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2.1 X-ray position of the nucleus 

M33 X-8 is severely affected by pile-up in the two 50 ks ACIS 
observations. The count rate in the ACIS-S observation is 0.5 
counts/frame compared to 3.0 counts/frame in the 10 ks non piled- 
up 10 ks ACIS-S observation (see §3 below). The pile-up fraction is 
therefore about 84% which makes it impossible to study the radial 
extension of the source to higher accuracy than the ROSAT HRI. 
Observations with the Chandra HRC would be able resolve any 
contribution from extended emission or faint point sources within 
5". The inner pixels of the ACIS-S (Fig. 1) are very piled-up hence 
have a dearth of counts so that the 0.5-2 keV image is doughnut 
shaped (standard G02346 events). Going to higher energies reduces 
the effect and there is a very well defined peak in both ACIS-S and 
I. 

CELLDETECT (CIAO v2.2.0.l|) locates M33 X- 
8 at (4098.1,4057.2)±0.1 pixels in the ACIS-I and at 
(4063.0,4026.3)±0.5 pixels in the ACIS-S. To check these 
values, we redefined the center as the location which minimizes 
the distances of all the events within a 300 pixel wide box roughly 
centered on M33 X-8 (this assumes circular symetry and avoids 
the problem posed by a pile-up doughnut). The values obtained by 
selecting only the events with energies above 1, 2 or 5 keV are all 
consistent with the CELLDETECT locations and errors which we 
therefore adopt. 



2.2 Optical position of the nucleus 

The optical position of the nucleus we adopt is the 2MASS po- 
sition which is tied to the Tycho reference frame. The (compact) 
nucleus is detected with J = 12.08 ± 0.04 in this survey. At 
a (J2000) l h 33 m 50:89 and S 30°39'36'/75 with an uncertainty 
of 0'.'15 this is, as far as we know, the best available absolute posi- 
tion. The position of the close 09III star was then found relative to 
the nu cleus using the UV H ST PC image described in D ubus et al. 
(1999b). We find a l n 33 m 50?85 and 5 30°39'37'.'65 to the 2MASS 
uncertainty. 

As an independent check, we cross-correlated the brightest 
stars in the various HST PC images of the nucleus ( Dubus et al. 
1999b) with the M33 star catalogue produced by the DIRECT 
project (Macri et al. 2001a). The DIRECT positions are tied to the 
USNO-A2 reference frame with an accuracy of / .'4-0 / /7. We found 
4-6 stars (depending upon the sky orientation of the PC chip) that 
correlated well with DIRECT positions, all lying at the same offset 
within 0'.' 1-0'.' 4 of each other. We found a mean corrected position 
for the nucleus of a l h 33 m 50:90 and S 30°39'37'.'l4. This is con- 
sistent with the 2MASS position within the limits of the DIRECT 
astrometric solution. 



2.3 X-ray source list and cross-correlation 

We applied CELLDETECT to identify X-ray sources in the two 
50 ks observations. The detection treshold was 4a, resulting in 5 1 
(77) possible sources in the ACIS-I (ACIS-S) dataset. We com- 
pared thi^hsrwidTthe^OSAT?^ & Pietsch 
(2001) to identify previously known X-ray sources. ROSAT and 
Chandra sources lying within 5" of each other were considered 
identical and are thereaf ter identified by their number in Table 1 of 
Haberl & Pietsch (200 1|). 
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Figure 1. Radial profiles of M33 X-8 in the two piled-up 50 ks ACIS-I 
(solid lines) and S (dashed line) observations (1 pixel=0.492"). All G02346 
events with energies > 0.5 keV are used in the top two profiles (9900 and 
6600 counts respectively for ACIS-I and S) while only those above 5 keV 
are used in the bottom profiles (3100 and 1200 counts; profiles are shifted 
down by a factor 10 for clarity). The profiles are (background subtracted) 
area-averaged number of events within x and x-1 pixels of the adopted cen- 
ters (§2.1). Pile-up appears most clearly as a flattening of the ACIS-S profile 
at low energies. 



We cross-co rrelated the X -ray source list wit h the Tycho-2 
( |Hog et al. 200C|), US NO-A2 ( |Monet et al. 1998| ) and 2MASS 
( JSkrutskie et al. 1997b catalogues. We also used a subset (M acri 



http://cxc.harvard.edu/ciao/ 



et al. 2001b; pVIochejska et al. 200l| ) of the DIRECT catalogue con 
taining variable opti cal sources (excludin g cepheids) and the radio 
source catalogue of Gordon et al. (1999). The latter contains opti- 
cally confirmed supernova remnants (SNR) in M33 and stellar-like 
background radio sources. The radio positions are accurate with re- 
spect to the ICRS frame to 7"/2<t, where 7" is the beam size and a 
is the signal-to-noise ratio of the source. 

We started by searching for all possible counterparts in the 
catalogues within 10" of each X-ray source and plotting their pixel 
offsets relative to the Chandra position. A given X-ray source can 
be associated with several sources from the same or different cata- 
logues. In the plot, true counterparts will cluster at the offset of the 
X-ray astrometry (the orientation error and ACIS chip distortion 
are negligible), while spurious identifications will appear uniformly 
distributed (the field of view is much greater than the searched cor- 
relation radius). The spread of the cluster is a measure of the ac- 
curacy of the astrometric solution. There was an obvious concen- 
tration of matches around the offset between X-8 and the 2MASS 
position of the nucleus (Table 2 and Fig. 2). Note that some of the 
matches may refer to the same object such as the USNO and Tycho 
match for HP87. None of the DIRECT variable sources matched an 
X-ray source within 3" although we note t hat a star with a 3.4 days 
optical period (D33J013334.8+30321 1.4, |Mochejska et al. 200l| ) 
lies within lfy_^Hhc_^ i 45_da^ i £dn3sn^ 
et al. 1999a). 
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Fi gure 2. Plot of the o ffsets in ACIS pixels (1 pixel=0.492' 

(o, ijordon et al. 1999) catalogues (see Table 1). The offset between the Chandra position of X 



The offset between X-8 and the UV bright star 
being associated with the nucleus (•) rather than with the nearby star (*) 



) between X-ray sources and their matches in the USNO (x), Tycho (+), 2MASS (o), and radio 

and the 2MASS position of the nucleus is indicated by a •. 
1" NNW of the nucleus is shown by a *. The radio counterparts cluster at an offset consistent with the ULX 



2.4 Position of the nuclear X-ray source 

The clustering of matches around the offset between the 2MASS 
position of the nucleus and the Chandra position of X-8 is partic- 
ularly clear for the radio correlations (Fig. 2). For uniformly dis- 
tributed sources, the probability p of having an X-ray source within 
d arcsec of a possible counterpart is p w (d/D) 2 if d <C D and 
where D is the radius of the field of regard. Correlating N X-ray 
sources with M possible counterparts in the field of view, the prob- 
ability P of observing C or more sources within d arcsec of a coun- 
terpart is 



tainty of C/.'6 (the mean weighted distance of the radio offsets for 
both ACIS-I and S). 

There is a radio source associated with the nucleus. This 
source (GDK102) has F 2 ocm = 0.6 ±0.1 mjy, F 6cm = 0.2 ± 
0.1 mjy is located C/.'6 away from the 2MASS position (and C//4 
away from the mean radio offset) and is a possible counterpart 
to M33 X-8. Using the results of the above analysis, we note the 
association of X-ray source HP84 with 2MASS 0133418+303849 
(J = 12.62 ± 0.02) is l ikely to be real. HP84 i s classified as a 
supersoft X-ray source by Haberl & Pietsch (2001). 



p = i-E(T)^ 1 -^ M "" < 

For C — 8 matches with the radio catalogue, d — 3" , N — 32, 
M = 126 and D « 8', the probability of seeing a cluster of 8 
associations within 3" is < 10 -11 . The same calculation for the 
ACIS-S observation (C = 7, N = 57, M = 117) yields a proba- 
bility <5 10~ 8 . We conclude the radio associations are very likely 
to be true counterparts (in contrast with the USNO, 2MASS and 
Tycho associations for which the same calculation does not rule 
out chance superpositions). 

We now estimate the probability that the offsets given by 
the radio counterparts are consistent with the offset between the 
2MASS position of the nucleus and X-8 (• in Fig. 2) or the nearby 
star position and X-8, plotted as a * in Fig. 2). A x 2 test on the 
measured distances between the radio and nucleus offsets yields 
a probability 0.28 (0.53) that these are consistent for the ACIS-I 
(ACIS-S) data. The probability for the radio offsets to be consistent 
with the nearby star are 2 • 10 -6 (ACIS-I) and 4 • 10~ 4 (ACIS-S). 
The ACIS-S result is less conclusive than the ACIS-I because of 
the higher uncertainty in the position of X-8 in the ACIS-S (±0.5 
pixel compared to ±0.1 pixel, see §2.1). We conclude the nearby 
star is ruled out as a possible counterpart at the 4.6a level and that 
M33 X-8 is at the 2MASS position of the nucleus with an uncer- 



3 X-RAY SPECTRUM 

We performed spectral analysis using data from the ACIS-S ob- 
servation Obsld 787. We extracted 27259 counts from a circle 
within 20" of the source position for a total of 3.06 counts/frame. 
The source was 7.7' off-axis, at which the the 50% enclosed en- 
ergy function (EEF) at 1.5 keV is 4" vs. 0.3"on-axis, resulting 
in the source covering an area larger by a factor ~170, to 0.02 
counts/frame/PSF, at which the estimated pileup fraction should 
be negligible (<1%). We confirmed this by comparing the distri- 
bution of photon event grades (for grades 0, 2, 3, 4 and 6) below 
2 keV and in the 2-8 keV range with those of a non-piled up X- 
ray observation (of ZW3146, Obsld 1651, Seq. Num 800119) and 
found that they were roughly consistent, though not identical; how- 
ever, both were similar to the branching ratio observed from a non 
piled-up observation of PG 1634-706 n, and so pileup does not ap- 
pear to be significant during this observation. Background counts 
were taken from an annulus of inner- and outer- radii of 25" and 
100". The expected number of background counts in the source re- 
gion is 167±3. The background subtracted source countrate is then 
2.92+/-0.02 c/s. 



3 http://asc.harvard.edu/cal/Hrma/hrma/misc/oac/psf2/index.html 
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Table 2. All possible matches within 3" of X-ray sources (Fig. 2). 



M33 X-8: Chandra ACIS-S/BI, Jan 1 1 2000 



Source 


Cat. 


a 


6 


AX 


AY 


err 








(degrees, J2000) 


(pixels) 


(") 




ACIS 


■I observation (Obsld 1730) 












HP102 


2MASS 


23.46205 


30.66021 


-0.0 


-2.4 


0.16 






GDK102 


23.46204 


30.66037 


0.0 


-1.2 


0.60 




HP84 


2MASS 


23.42445 


30.64703 


-0.1 


-2.0 


0.16 


'a 


HP74 


GDK64 


23.39958 


30.60785 


0.7 


-1.5 


0.26 


> 
o 


HP130 


GDK148 


23.54462 


30.70671 


-1.2 


-2.1 


0.40 




HP62 


USNO 


23.37106 


30.70483 


-2.0 


-1.9 


0.19 






GDK50 


23.37096 


30.70482 


-1.4 


-1.9 


0.42 




CXOl 


2MASS 


23.34009 


30.65008 


1.5 


0.1 


0.46 






2MASS 


23.34066 


30.64976 


-2.1 


-2.3 


0.46 






GDK38 


23.34046 


30.65010 


-0.8 


0.2 


1.26 




HP 106 


GDK112 


23.47858 


30.55301 


-0.2 


-1.9 


0.26 




HP80 


GDK67 


23.40621 


30.78870 


0.6 


-3.2 


0.52 




HP67 


GDK57 


23.38008 


30.55939 


-0.9 


-1.3 


0.27 




ACIS 


-S observation 


(Obsld 786) 












HP102 


2MASS 


23.46205 


30.66021 


0.6 


-1.2 


0.29 






GDK102 


23.46204 


30.66037 


0.6 


-0.0 


0.65 




HP98 


GDK100 


23.46025 


30.63937 


1.3 


-1.2 


0.25 




HP84 


2MASS 


23.42445 


30.64703 


0.8 


-0.8 


0.15 




HP62 


GDK50 


23.37096 


30.70482 


-0.1 


-1.4 


0.42 






USNO 


23.37106 


30.70483 


-0.8 


-1.3 


0.20 




HP106 


GDK112 


23.47858 


30.55301 


0.4 


-0.0 


0.22 




HP87 


Tycho 


23.43076 


30.77521 


1.5 


-0.3 


0.36 






2MASS 


23.43076 


30.77527 


1.5 


0.1 


0.38 






USNO 


23.43081 


30.77521 


1.2 


-0.3 


0.38 






USNO 


23.43013 


30.77541 


5.5 


1.2 


0.38 




CX02 


GDK51 


23.37192 


30.76057 


0.4 


-3.0 


0.41 




HP80 


GDK67 


23.40621 


30.78870 


1.2 


0.0 


0.48 




CX03 


2MASS 


23.49518 


30.83467 


-0.5 


2.8 


0.51 




HP64 


GDK52 


23.37279 


30.81954 


-0.5 


1.0 


0.94 





HP refers to the Haberl & Pietsch (2001) catalog of ROSAT sources and 



is used to identify the sources (column 1). CXO refers to Chandra sources 
which were not previously detected. GDK refers to the jordon et al. (199S ) 



catalog of radio sources in M33 (their Table 1). For each match we list the 
catalogue (column 2) and position in this catalogue (columns 3 & 4). AX 
and AY (columns 5 & 6) are the offsets in ACIS pixels between the X-ray 
and matched source position (plotted in Fig. 2). Column 5 is the offset error 
(cr 2 , + o" 2 at ) 1//2 , where cr 2 , is the error in the X-ray source position (given 
by CELLDETECT, see §2.3) and o- 2 at is the error in the position of the 
correlated source (given in the relevant catalogue, see §2.3). 



Using XSPEC v 11.0.1 ( [Arnaud 199rj| ), we produced a pulse- 
height-analyser (PHA) spectrum of the extracted counts data be- 
tween 0.6-9 keV, binning at 130 eV between 0.6 and 2.0 keV, at 
260 eV between 2.0 and 5.4 keV, and with two energy bins be- 
tween 5.4 and 9.0 keV. We included a 5% systematic uncertainty in 
the fits. Best-fit parameters for acceptable models are in Tableg. 

An absorbed powerlaw spectrum (wabs *powerlaw) is re- 
jected, with a value of reduced chi-squared statistic x^=4.12 for 
23 degrees of freedom (dof), with a null-hypothesis probability of 
10 -10 . The power-law spectrum over-predicts counts in the 2-4 
keV range, and under-predicts them in the 4-9 keV range. An ab- 
sorbed blackbody spectrum is also statistically rejected (xt = 21.4 
for 23 dof). 

We find a moderately acceptable fit (xS=l-35) for an ab- 
sorbed thermal bremsstrahlung spectrum. An absorbed disc black- 
body spectrum produces a better fit to the data (xS=0.95), with an 




Energy (keV) 

Figure 3. Best-fit absorbed disc-blackbody model (solid-line histogram), 
with the observed Chandra/ ACIS-S spectrum of M33 X-8 (crosses), cor- 
rected for the detector response, in vF v . The deviation of the data points 
near 3 keV from the model are not significant (2<r); the excess near 6-9 keV 
deviates from the best-fit model at the 3<r level, and is consistent with pos- 
sible contributions from pile-up at the ~ 1 % level, which are not accounted 
for in the best-fit spectrum. 



Table 3. Chandra Observation Spectral Parameters (0.6-9 keV) 



Thermal Bremsstrahlung 

N H (10 22 cm~ 2 ) 
kT 

J n e njdV 
xl (dof) 

Disc Blackbody 

N H (10 22 cm- 2 ) 

71 
in 

RinV cos 6 
xl (dof) 



0.19±0.01 
3.5±0.1 keV 
(1.74±0.05) xlO 62 cm- 
1.35 (23) 



().06±0.01 
1.18±0.02keV 
57±3 km (D/795 kpc) 
0.95 (23) 



Assumed source distance d = 795 kpc. Uncertainties are la. 

X-ray column density Njj — (6± 1) x 10 20 cm -2 , consistent with 
galactic absorption only Q The disc blackbody inner temperature 
fcTi n and inner-radius RinV cos 9 ar e identical to the va lues found 
in earli er observations with ASCA (Takano et al. 1994) and Bep- 



poSAX (Parmar et al. 2001). The flux (corrected for absorption) of 
the best-fit spectrum is 2. Ox 10" 11 erg cm -2 s _1 , for a luminos- 
ity 1.5xl0 39 erg s _1 (0.5-10 keV), consistent (including spec- 
tral uncertainty) with previously measured values. Unlike previous 
work, a power-law component is not required by the data. When we 
include a power-law component in addition to the absorbed disc- 
blackbody model, holding a — 2 fixed, the 90% confidence upper- 
limit for the flux of this component is 1.1 x 10 _11 erg cm -2 s _1 . 



4 X-RAY VARIABILITY 

We examined the data for long-term noise-like time variability. We 
binned the 27259 counts detected in the full Chandra pass band, 
with 1.04104 sec time resolution, and took a Fourier transform 



W3NH at ittp://heasarc.gsfc. nasa.gov/ (Dickey & Lockman 1990) 
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of the data. With the resulting Four ier components, w e obtained 



a power-density spectrum (PDS; see Lewin et al. 1988 for the use 
of PDS for broad-band timing analysis) which we rebinned loga- 
rithmically. Examination of the resulting PDS shows no evidence 
of power above that expected from Poisson noise. Modelling the 
power as a power-law of slope=l plus a constant (held fixed at the 
Poisson level) in frequency (P(i>) = A v~ 1 + C), we find an upper- 
limit to the root-mean-square variability of <3% (90% confidence; 
0.0001-0.5 Hz). 

When we binned the data by a factor of 1, 10 and 100, (in 
which the average counts per bin are 2.938, 29.4, and 294 counts), 
we find that the maximum counts per bin was 11, 48 and 330 
counts, which is consistent with the maximum number of counts 
expected from Poisson fluctuations. Thus, we find no evidence for 
intensity variability on timescales between 1 and 10,000 sec. 



5 DISCUSSION 

The high luminosity of M33 X-8 is not sufficient to rule out a neu- 
tron star X-ray binary since magnetospheric accretion can result in 
super-Eddington luminosities. This was observed e.g. in the high 
mass X-ray binaries (HMXBs) SMC X-l or A0535-66. A HMXB 
would have been an attractive possibility should M33 X-8 have 
been associated with the nearby 09III star. The Chandra localisa- 
tion of X-8 rules this out (§2). X-8 could still be a HMXB hiding in 
the nucleus but this would be surprising since the compactness and 
spectra of the nucleus make it u nlikely that it hosts > IOMq stars 
( Long, Charles & Dubus 2002| ); a HMXB evolves on the uncom- 
fortably short thermal timescale of the donor star (< 10° yrs for 
M2 > IOA'/q); the X-ray spectrum of X-8 is unlike that of typical 
HMXBs which have flat power laws in the 0.5-10 keV range; and 
the power spectrum shows no evidence for pulsations. 

On the other hand, the X-ray spectrum and absence of short- 
term X-ray variability below 1 Hz are highly reminiscent of the 
high or very higj^s^ate^rf^alacticjs^ 

6 Lewin 1995). In both states the spectrum is dominated by a disc 
blackbody with a weak contribution from a steep power law. The 
very high state shows 3-10 Hz QPOs with a few % rms superposed 
on 1-10% rms flat top noise in the 1-10 Hz range. There is little 
variability in the high state. GRS 1915+105, a microquasar in our 
Galaxy w ith a ~ 14 Mq black hole accreting f rom a ~ 1.2 Mq K 
giant star ( Greiner, Cuby & McCaughrean 2001 ), also has such soft, 
low variability states ( Belloni et al. 2000| ). The optical emission 
from such a system would be lost in the nucleus of M33. Since its 
discovery in 1992, GRS 1915+105 has been steadily detected with 
X-ray luminosities around or in excess of 10 39 erg s _1 . 

The analogy with GRS 1915+105 is even more attractive 
considering the 0.6 mjy (respectiv ely 0.2 mJy) radio source de- 
tected by the VLA at 20 cm (6 cm; Gordon et al. 1995 ) which we 



show to be consistent with the ultra-luminous X-ray source (§2.4). 
GRS 1915+105 has been known to show steady radio emission 
around 30-50 mjy at 2 cm, framed by powerful ejection episodes 
where the emission from the relativistic jets becomes > 500 mjy 
(Fender et al. 1999). The radio source in M33 is more powerful 
than that; but a small change in the beaming angle of the jet (~ 70° 
from the line-of-sight for GRS 1915+105) could easily increase the 
observed radio flux. Follow-up observations of the radio source are 
needed to test this hypothesis. 
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